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Abstract 

We perform Hartree calculations of symmetric and asymmetric semi-infinite nuclear 
matter in the framework of relativistic models based on effective hadronic field theories 
as recently proposed in the literature. In addition to the conventional cubic and 
quartic scalar self-interactions, the extended models incorporate a quartic vector self- 
interaction, scalar-vector non-linearities and tensor couplings of the vector mesons. 
We investigate the implications of these terms on nuclear surface properties such as 
the surface energy coefficient, surface thickness, surface stiffness coefficient, neutron 
skin thickness and the spin-orbit force. 
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1 Introduction 



Quantum hadrodynamics (QHD) and the relativistic treatment of nuclear systems has 
been a subject of growing interest during recent years The a— u model of Walecka 

[|l| and its non-linear extensions with cubic and quartic self-interactions of the scalar-meson 
field have been widely used to this end. This model contains Dirac nucleons together 
with neutral scalar and vector mesons as well as isovector-vector p mesons. At the mean 
field (Hartree) level, it already includes the spin-orbit force, the finite range and the density 
dependence which are essential ingredients of the nuclear interaction. This simple model has 
become very popular in relativistic calculations and describes successfully many properties 
of the atomic nucleus. 

From a theoretical point of view, the non-linear o — uj model with cubic and quartic 
scalar self-interactions was classed within renormalizable field theories which can be char- 
acterized by a finite number of coupling constants. However, very recently, generalizations 
of this model that include other non-linear interactions among the meson fields and tensor 
couplings have been presented on the basis of effective field theories by Serot et al. 
The effective theory contains many couplings of non-renormalizable form that are consis- 
tent with the underlying symmetries of QCD. Consequently, one must find some suitable 
expansion parameters and develop a systematic truncation scheme. For this purpose the 
concept of naturalness has been employed: it means that the unknown couplings of the 
theory should all be of the order of unity when written in appropriate dimensionless form 
using naive dimensional analysis H^HH • Then, one can estimate the contributions coming 
from different terms by counting powers in the expansion parameters and truncating the 
Lagrangian at a given level of accuracy. 

One important fact is the observation that at normal nuclear densities the scalar and 
vector meson fields, denoted by <£> and W, are small as compared with the nucleon mass 
M and that they change slowly in finite nuclei. This implies that the ratios $/M, W/M, 
|V$|/M 2 and |VW|/M 2 are useful expansion parameters when the effective field theory is 
applied to the nuclear many-body problem. From this viewpoint, if all the terms involving 
scalar and meson self-interactions are retained in the Lagrangian up to fourth order, one 
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recovers the well-known non-linear a — uj model plus some additional terms For the 

truncation to be consistent, the corresponding coupling constants should exhibit naturalness 
and cannot be arbitrarily dropped out without an additional symmetry argument. The 
effective Lagrangian truncated at fourth order contains thirteen free parameters that have 
been fitted to reproduce twenty-nine finite nuclei observables |5]||. Remarkably, the fitted 
parameters turn out to be natural and the results are not dominated by the last terms 
retained. This evidence confirms the utility of the principles of naive dimensional analysis 
and naturalness and shows that truncating the effective Lagrangian at the first lower orders 
is justified. 

The term with a vector-meson quartic self-interaction has been considered previously in 
relativistic mean field (RMF) calculations from a phenomenological point of view. Bodmer 



10 1 considered this coupling to avoid the negative coefficient of the quartic scalar self- 



interaction that appears in many non-linear a-w parametrizations that correctly describe the 
atomic nucleus [0]] . In some special situations this negative term can lead to a pathological 
behaviour of the scalar potential. On the other hand, the equation of state is softened at 
moderate high densities when the vector non-linearity is taken into account. The quartic 



vector self-interaction has also been phenomenologically used by Gmuca [fT^Jl3f| in a non- 
linear a—uj model for parametrizing Dirac-Brueckner-Hartree-Fock calculations of nuclear 
matter. The same idea was developed by Toki et al. and applied to study finite nuclei 



14| and neutron stars ||15|| . Recently, the properties of high-density nuclear and neutron 



matter have been analyzed in the RMF approach taking into account scalar and vector 
non-linearities ||. 

The tensor couplings of the vector uj and p mesons to the nucleon were investigated by 
Reinhard et al. |3],[16| as an extension of the RMF model, and more recently by Furnstahl 
et al. P|i7[| from the point of view of relativistic effective field theory. In these works it was 
shown that the tensor coupling of the u meson has an important bearing on the nuclear 
spin-orbit splitting. 

The surface properties of nuclei play a crucial role in certain situations. This is the 
case, for instance, of saddle-point configurations in nuclear fission or fragment distributions 
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in heavy-ion collisions. Within a context related to the liquid droplet model (LDM) and 



the leptodermous expansion [fig] , the surface properties can be extracted from semi-infinite 
nuclear matter calculations either quantally or semiclassically (though the total curvature 
energy coefficient can only be computed semiclassically Jim). In the non-relativistic case 
most of the calculations of the surface properties have been carried out using Skyrme forces, 
quantally |2(J or semiclassically with the help of the extended Thomas-Fermi (ETF) method 
PU| , PT| . In the relativistic case the nuclear surface has been analyzed within the a— u model 
since a long time ago. The calculations have been performed semiclassically using the rela- 
tivistic Thomas-Fermi (TF) method or its extensions (RETF), for symmetric fl|,|6|,f22| p3] , p4T 
and asymmetric [ P5|J2B[ ] matter, and also in the quantal Hartree approach [ r27|J25|J2P| ] . 

In the framework of the relativistic model and effective field theory, the main purpose 
of the present work is to carefully analyze the influence on surface properties of the quartic 
vector non-linearity, of the newly proposed scalar-vector self-interactions and of the tensor 
coupling. We shall investigate quantities such as the surface energy coefficient, surface thick- 
ness, spin-orbit strength, surface stiffness coefficient and neutron skin thickness obtained 
from Hartree calculations of symmetric and asymmetric semi-infinite nuclear matter. 

The paper is organized as follows. Section 2 is devoted to the basic theory. The results 
on the surface properties of symmetric matter are discussed in Section 3. Section 4 addresses 
the case of asymmetric systems. The summary and conclusions are given in the last section. 
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2 Mean field equations for symmetric semi-infinite nu- 
clear matter 



Following Ref . , to derive the mean field equations one starts from an energy functional 
containing Dirac baryons and classical scalar and vector mesons. The energy functional 
can be obtained from the effective Lagrangian in the Hartree approach using many-body 
techniques |5]||. However, this energy functional can also be considered as an expansion in 
$/M, W/M, |V$|/M 2 and |VW|/M 2 of a general energy density functional that contains 
all the correlation effects. The theoretical basis of this functional lies on the extension of 
the Hohenberg-Kohn theorem |30| to QHD [|3T|. Using the Kohn-Sham scheme [3^| with 



the mean fields playing the role of Kohn-Sham potentials, one finds similar mean field 
equations to those obtained from the Lagrangian f3lH , but including effects beyond the 
Hartree approach through the non-linear couplings 



A semi-infinite system of uncharged nucleons corresponds to a one-dimensional geometry 
where half the space is filled with nuclear matter at saturation and the other half is empty, so 
that a surface develops around the interface. The fields and densities change only along the 
direction perpendicular to the medium. Specifying the energy density functional considered 
in Refs. || and |J to symmetric semi-infinite nuclear matter with the surface normal pointing 
into the z direction one has 

= E^)^-V + m-^)]+W(z)-^/3oc-VW(z)^ a (z) 

A 1 + ai— — (V$(zi) + - + — — — + — — — $ (z) 

2g 2 B [ L ^ ai m J { { >> \2 3! M 4! M 2 J g y K ' 

-i(l + ^ + f^)f^). (-) 

where the index a runs over all occupied states of the positive energy spectrum, $ = g s (p 
and W = g v Vo (notation as in Ref. JIJ). Except for the terms with ot\ and «2, the functional 
(|2.1| ) is of fourth order in the expansion. We retain the fifth-order terms ot\ and «2 because 



in Refs. || and they have been estimated to be numerically of the same magnitude as 
the quartic scalar term in the nuclear surface energy. 

The mean field equations are obtained by minimizing with respect to (p' a , $ and W: 

i-ia- V + (3[M - $(*)] + W(z) - % Jz.f3oL-VW{z)\ <p a (z) = e a <p a (z) , (2.2) 



m l^2 ( ..\ ( K 3 , ^ Hz) 



M w \ 2 3! M J 
2M\ n 12 M J K ' 

- ^CoW 3 (z) + g[VS(*) • W(z) + $(z)A^)] . (2.4) 
The baryon, scalar and tensor densities are respectively 

M = 5>i(*)v«(*) , (2-5) 

= I>£(*)/W*), (2-6) 

M*) = EtjV-[^W«^)]- (2.7) 

a 

The expression of the four-component spinors (p a (z) in the semi-infinite medium was given 
by Hofer and Stocker in Ref. [ pTj] , 

In a semi-infinite nuclear matter calculation the sum over the single-particle states is 
replaced by an integration over momenta: 

?^?/«*. (-> 

where Q stands for the volume of the box, the factor 2 takes into account the isospin degree 
of freedom and A describes the spin orientation of the nucleons. Introducing the Fermi 
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momentum hp, the integration domain is restricted to k 2 . + k 2 + k 2 = k\ + k 2 < kp, with 



k z > if the bulk nuclear matter is located at z = — oo. Following the method outlined in 
Ref. [27 1 one finds two sets (A = ±1) of first-order differential equations for the orbital part 
of the upper and lower components of the Dirac spinors: 



dGgjz) 
dz 

dF a (z) 
dz 



Xk± + 



Xk ± + 



fy dW{z) 

2M dz 

fy dW(z) 

2M dz 



GJz) 



Fjz) 



[e a -W(z) + M*(z)]F a {z) 



(2.9) 



[e a -W(z)-M*(z)]G a (z), (2.10) 



where a = (k z ,k±, A) and M*(z) = M — &(z) is the Dirac effective mass of the nucleons. 
From the asymptotic behaviour at z = — oo (bulk nuclear matter), the condition on the 



energy eigenvalues is e a = y k\ + k 2 + M^ 2 + Woo, with and Woo being the nuclear 
matter values of M* and W. The densities for each spin orientation A = ±1 read 

2 r k F 



7T 2 JO 



dk 7 



2 f k v 

7T Z JO 



V^F k z 



dk ± k ± (\G a (z)\ 2 +\F a (z)\ 2 ) 
dk ± k ± (\G a (z)\ 2 -\F a (z) 



/4(*) = — / dk s r dk ± k ± -(—F a (z)G a (z] 
ti z Jo Jo dz \M 



and the total densities are given by 



p(z) = J2p X (z) 



p 8 (z)=Y,Ps(z), PtW = £pt( 



(2.11) 
(2.12) 
(2.13) 

(2.14) 



Using the equations of motion the energy density of the semi-infinite nuclear matter 
system can be written as follows: 



S(z) 
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Finally, the surface energy coefficient E s is obtained from the expression [18 



/oo 
dz [E{z) - (ov + M)p{z)} 
-oo 



(2.16) 



where a v is the energy per particle in bulk nuclear matter and vq is the nuclear radius 
constant: r = (3/47rpo) 1 ^ 3 ; with p the nuclear matter density. 

Another important quantity in the study of the nuclear surface structure is the spin-orbit 
interaction. By elimination of the lower spinor in terms of the upper spinor, one obtains a 
Schrodinger-type equation with a term V so (z) that has the structure of the single-particle 



spin-orbit potential for the non-relativistic case |27|j29| . In our present model the orbital 
part of the spin-orbit potential reads as 



V so (z) 



1 



2M 



'dW(z) d$(z) 



+ 



/v dW(z) 



(2.17) 



e a - W(z) + M*(z) \ dz ' dz J ' M dz 
In the non-relativistic limit, by means of a Foldy-Wouthuysen reduction, Eq. ( [2.170 becomes 

1 



AM 2 



dz 



dz 



(2.18) 



and the nucleons are then moving in a central potential of the form 



V c {z) = W(z) - $(z) . 



(2.19) 
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3 Surface properties in the symmetric case 



Although previous works |^J^J^ , |^J27| , P8| , |29|| have thoroughly investigated the prop- 



erties of the nuclear surface in the standard non-linear a — to model, for which k 3 and k 4 
are the only non-linearities in fl2.1|) , here we wish to enlarge this study by including the 
additional non-linear and tensor couplings considered in Refs. 0,0,0. In concrete, we want 
to study the role of the quartic vector self-interaction ( that has been used after the work 
of Bodmer [|1IJ, the role of the terms with rji and 772 that couple the scalar and vector 
fields, of the terms with a\ and a 2 that imply the gradients of the fields, and of the tensor 
coupling f v of the vector u> meson to the nucleon. While ( , rji and r] 2 can be classified as 
volume contributions, the couplings a±, a 2 and f v are genuine surface terms. On the basis 
of the concept of naturalness there is no reason to omit any of these terms in the energy 
density functional ( |2.1| ), unless there exists a symmetry principle to forbid it. However, to 
clarify the impact on the surface properties of the aforementioned couplings we will analyze 
each one separately, as it has been similarly done in Ref. ]?]. In this section we shall study 
symmetric systems, while in Section 4 we shall address the case of asymmetric matter. 



3.1 Effect of the quartic vector self- interaction 

In the conventional non-linear a — uo model the value of the coefficients g^/m^, g^/m^, 
n 3 and K4 can be univocally obtained by imposing that for nuclear matter at saturation the 
density po, energy per particle a v , effective mass M^/M and incompressibility modulus K 
take given values. When the vector-meson quartic self-interaction is switched on, the Dirac 
equation for the baryons and the Klein-Gordon equation for the vector field in infinite 
nuclear matter become (throughout this subsection we set rji = 772 = «i = a 2 = f v = 0): 

Ov = y/kl + M^ + W^-M, (3.1) 

rr^Woo = S>o-~CoV^. (3.2) 
o 

The saturation density po and the Fermi momentum hp are related as usual by po — 2/cp/37r 2 . 
Specifying po, a v and M^/M, from the above equations one extracts the coupling constant 
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g v as a function of Co (the nucleon and u masses take their empirical values: M = 939 MeV 
and m v = 783 MeV). The steps to calculate g^/rn^, K3 and k 4 are then the same as when 



Co = 0, see e.g. Refs. and |IT| , but now these coefficients become functions of Co- The 
reader will find a detailed study of the implications of the vector non-linearity Co i n nuclear 
matter in Refs. |7j and |L0| . 

The assumption of naturalness requires that the couplings g s /47r, g v /47r, k 3 , k 4 and 
Co should all be roughly of the order of unity. Figure 1 illustrates the variation of these 
couplings as a function of the non-dimensional parameter 



used by BodmerQ in Ref. [TIJ. With m s = 490 MeV, the figure presents the results for four 
sets of nuclear matter properties: po — 0.152 fm~ 3 (kp = 1.31 fm -1 ), a v = —16.42 MeV, 
K = 200 and 350 MeV, and M*JM = 0.6 and 0.7. When K = 200 MeV and M^/M = 0.6 
the equilibrium properties of the interaction and the scalar mass m s are very close to those 
of the non-linear parametrization NL1 |33] . 



We realize that g s and g v are only weakly affected by r] (the change is not appreciable 
in the scale of Figure 1). However, r/o has a direct effect on the scalar-meson quartic self- 
interaction K4, moving it from a negative value when the quartic vector term is absent 
(7/0 = 00) to a desirable positive value when i] ~ 2. The change of sign of k 4 takes place 
at larger values of r] for larger K and M^/M; in fact, /t 4 is already positive at rj = 00 if 
K = 350 MeV and M^/M = 0.7. Both non-linearities of the scalar field k 3 and k 4 remain 
in the natural zone for 2 < i]q < 00 approximately, excepting the case of K = 350 MeV 
and M^/M = 0.7, but they start to depart appreciably from their natural values when 
7/0 < 2. These trends fairly agree with the assumption of naturalness: t]q > 2 corresponds 
to the region where Co can be considered as natural (see the figure), and for this range of t]q 
values also the rest of coupling constants are natural. The behaviour gleaned from Figure 
1 is rather independent of the saturation density p and energy a v , and e.g. we have found 



similar trends with the specific nuclear matter properties used by Bodmer in Ref. [10 



L Notice that in Ref. p0| the parameter 770 was called z. 
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Next we turn our attention to the surface properties. In Figures 2 and 3 we have plotted, 
respectively, the surface energy coefficient E s and the surface thickness t of the semi-infinite 
density profile (the 90%-10% fall-off distance) against the tjq parameter. The selected values 
of ?7o are those employed in Table 1 of Ref. With p = 0.152 fm -3 and a v = —16.42 
MeV, we have performed the calculations for a few incompressibilities (K = 125, 200 and 
350 MeV) and effective masses (M^/M = 0.6 and 0.7). Furthermore, we have considered 
two values of the mass of the scalar meson (m s = 490 and 525 MeV). This quantity governs 
the range of the attractive interaction and determines the surface fall-off: a larger m s results 
in a steeper surface and a reduction of E s and t. 

From Figures 2 and 3 we realize that the quartic vector self-interaction scarcely alters the 
values of the surface energy coefficient and of the surface thickness if remains in the natural 
domain, i.e., if 2 < 770 < 00. At K — 350 MeV and M*/M > 0.6, E s is raised by decreasing 
770- This tendency may be inverted at smaller values of the incompressibility, depending also 
on the value of the nucleon effective mass, but the global trends are practically independent 
of the mass of the scalar meson. The thickness t exhibits a more monotonous behaviour: 
in all cases it stays almost equal to its value at 770 = 00 and goes down slightly for 770 < 2. 
Noticeable departures of the surface energy coefficient and thickness from the r] = 00 
values can be found only if 770 is decreased beyond the natural limit, a situation where the 
interaction is mainly ruled by the vector-meson quartic self-interaction. 

Even though the impact of the vector-meson quartic self-interaction is small, it can help 
to find parameter sets for which both the surface energy coefficient and the surface thickness 
lie in the empirical region. This fact is illustrated in Figure 4, where E s and t are drawn 
versus M^/M for 770 = 00, 2 and 0.5 and for m s = 450, 500 and 550 MeV, using the nuclear 



matter conditions of Ref. p = 0.1484 fm~ 3 , a v = -15.75 MeV and K = 200 MeV. 
(We also performed the calculations for 770 = 5 and 770 = 1 to confirm the trends we discuss 
below.) The horizontal dashed lines in Figure 4 serve to indicate the empirical region for 
the surface energy and thickness. 

If at 770 = 00 we concentrate, for instance, on the parametrizations with m s = 450 MeV 
we see that the one with M^/M = 0.7 yields, simultaneously, E s and t within the empirical 
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region. The agreement between the calculated values with m s = 450 MeV and the empirical 
region improves for rj = 2, where both E s and t are acceptable for practically all the values 
of the effective mass considered. 

For 1] > 2 the dependence of E s and t upon the nucleon effective mass is similar to that 



found in the usual a— u model without a quartic vector self- interaction |23| , P8| . The general 
tendencies start to change when t]q is lowered and leaves the natural region. As rjo becomes 
smaller Figure 4 shows that the slope of the curves of E s and t as a function of M^/M 
changes, and that the curves for the different m s come closer together. 

To get more insight about the influence of the vector-meson quartic self-interaction, we 
display in Figure 5 the profiles of the baryon density p(z) and of the surface tension density 
(Swiatecki integrand) a(z) = £(z) — (a v + M)p(z), Eq. Q2.16|) . In turn, we have represented 
in Figure 6 the orbital part of the spin-orbit potential V so (z), Eq. (|2.17|) , at the Fermi 
surface (i.e., evaluated at k = kp) and the central mean field V c (z) defined in Eq. ( |2.19p . 
The properties of the interactions used in these figures are po = 0.152 fm -3 , a v = —16.42 
MeV, K = 200 MeV, M^/M = 0.6 and 0.7, and m s = 490 MeV. Results are shown for 
rjo = oo, 2 and 0.5. The corresponding values of the surface energy coefficient and surface 
thickness can be read from Figures 2 and 3. 

It can be seen that the local quantities depicted in Figures 5 and 6 oscillate as functions 
of z (Friedel oscillations). Both the surface tension density and the single-particle spin- 
orbit potential are confined to the surface and average to zero inwards as the bulk matter is 
approached. The local profiles for r] = 2, which somehow marks the limit of naturalness as 
we have commented, are almost equal to those obtained in the absence of the quartic vector 
self- interaction. In agreement with Figures 2 and 3, only when rjo is decreased to non-natural 
values one can notice some changes in the profiles, which are more visible for M^/M = 0.6 
than for M^/M = 0.7. Decreasing r] Q makes the surface steeper and the thickness t smaller, 
produces an enhancement in the surface region of the density p(z) and of the mean field 
V c (z), and builds up Friedel oscillations in a(z) and in the spin-orbit potential V so (z). 

It is well known that the experimental spin-orbit splittings require within narrow bounds 
a Dirac effective mass M^/M around 0.6 in the conventional relativistic model |TT|j7]1. Figure 
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6 shows that introducing a quartic vector self-interaction makes the spin-orbit well deeper. 
However, it is only a minor effect: with M^/M = 0.7 it is not possible to reproduce the spin- 
orbit interaction of the case M^/M = 0.6 at 770 = 00, not even if one sets 770 = 0.5 (which 
in addition brings about an unreallistically small t, see Figure 3). A similar conclusion 
was drawn in Ref. [10[ from an analysis in nuclear matter. We also have computed the 
non-relativistic limit V™(z) of the spin-orbit potential given by the expression ( |2.18 ). 
In agreement with Ref. we have found that while V™(z) qualitatively reproduces the 
behaviour of V so (z), it strongly underestimates the quantitative depth of the fully relativistic 
spin-orbit strength (by ~ 20% for M^/M = 0.6 and by ~ 15% for M^/M = 0.7). 

3.2 Influence of the volume cubic and quartic scalar-vector inter- 



The next bulk terms in the energy density ( |2.1| ) that contain non-linear meson interac- 
tions are 



In analogy to Figure 1 for r/ Q , in Figure 7 we study the change of the couplings g s /4:7i, g v /47r, 
K3 and K4 with the parameters rji and 772 (introducing each one separately) for some specific 
equilibrium properties. With p = 0.152 fm -3 , a v = —16.42 and m s = 490 MeV, in part (a) 
of Figure 7 it is K = 200 MeV and M*JM = 0.6, while in part (b) it is K = 350 MeV and 
M*JM = 0.7. 

For K = 200 MeV and M^/M = 0.6 both rji and 772 have a considerable effect on the 
scalar non-linearities k 3 and k 4 . The coupling k 4 changes sign for 772 ~ 1 , but it remains 
negative in the interval of r)i values used. To keep all the coupling constants within natural 
values we see that the range for 771 and 7/2 (when introduced separately) is restricted to 
run roughly from —0.5 to 2.5. Consider now different values of the incompressibility and 
effective mass, as in part (b) of Figure 7. The dependence of the couplings on 772 is not 
significantly altered. However, increasing either K or M^/M results in a smoother slope of 
K3 with 771, while it makes K4 grow steadily with rji and become positive at some value of 
this parameter. We have checked for K = 200 MeV and M^/M = 0.55, and for K = 125 



actions 




(3.4) 
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MeV and M^/M = 0.6, that k 4 remains negative at all values of 771 and that k% grows with 
r)i faster than in part (a) of Figure 7. 

The influence of the non-linear interactions rji and 7] 2 on the surface properties is analyzed 
in Figure 8. With ( = a.\ = a 2 = f v = 0, in this figure we have computed E s and t taking 



into account the terms (|3.4j) . The saturation conditions of the interaction and the scalar 
mass are the same as in part (a) of Figure 7. The results are displayed in the plane 771-772 in 
the form of contour plots of constant E s (solid lines) and of constant t (dashed lines). The 
range of variation of rji and 772 lies in the region imposed by naturalness, and yields values 
of E s and t within reasonable limits. 

As it can be inferred from the nearly vertical lines in the 771-772 plane, the surface energy 
coefficient and thickness depend mostly on 771 and are rather independent of 772. The con- 
sequence of increasing 771 is a reduction of the values of E s and t. The lines of constant t 
turn out to be, roughly speaking, parallel to the lines of constant E s . This means that from 
the interplay of the parameters 771 and 772 it is not possible to change the value of t relative 
to that of E s (for example, we see in Figure 8 that t ~ 2.2 fm if E s = 18 MeV). We have 
calculated the spin-orbit potential V so (z) at the Fermi surface for several values of 771 and 
772. We have found that these couplings have a marginal effect on the spin-orbit strength, 
as it happened to be the case with the other bulk non-linearity 770. 

To get some information about the incidence on the surface energy and thickness of all 
the volume non-linear meson interactions together, we have repeated the calculations in the 
771-772 plane setting 7/0 = 2 for the quartic vector self-interaction. One finds similar features 
to those of Figure 8. The effect of 770 = 2 is just shifting E s and t towards smaller values as 
compared with the case 770 = 00 (Co = 0), which is in accordance with what was found in 
Figures 2 and 3 at K = 200 MeV and M^/M = 0.6. 

3.3 Influence of the non-linear terms with gradients 

Now we discuss the non-linear interactions 



1 $ 

2M 



% (V$) 2 - % (W) 2 



(3.5) 
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that vanish in infinite nuclear matter. We recall that these terms are actually of order 5 in 
the expansion of the effective Lagrangian but, following Refs. and |J, we include them 
because they can be relevant in the surface due to their gradient structure. 

Using the same saturation properties and scalar mass of Figure 8, in Figure 9 we have 
calculated E s and t for several values of a\ and ot 2 with (o = rji = r] 2 = f v = 0. One observes 
that the curves of constant E s are projected onto the plane a%-a 2 as almost parallel straight 
lines (at least in the analyzed region, corresponding to natural values of ol\ and a 2 \ The 
same happens to the curves of constant t. But in contrast with the situation found in the 
plane 771^72 (Figure 8), the slope of the lines of constant t is different from that of the lines 
of constant E s . This means that by varying a± and a 2 one can achieve some modification 
on the surface thickness while keeping the same surface energy. For example, if we consider 
the contour line of E s = 18 MeV we find that for a 2 = 2.0 it is t ~ 2.05 fm, whereas for 
a 2 = —1-5 it is t ~ 2.25 fm. From Figure 9 we also see that increasing «i at constant a 2 
brings about larger values of E s and t, and that the opposite happens if one increases a 2 at 
constant ot\. 

We have repeated the calculations of Figure 9 {K = 200 MeV, M^/M = 0.6) for 
K = 350 MeV and for M^/M = 0.7, to verify to which extent the behaviour in the ai-a 2 
plane is affected by the incompressibility and effective mass of the interaction. Certainly, 
the contour lines of E s and t are shifted with respect to Figure 9, but the trends with ai 
and a 2 turn out to be qualitatively the same. The range of variation of the surface energy 
and thickness in the a\— a 2 region we are considering is shorter when M^/M = 0.7, while 
it is more or less the same when K = 350 MeV. 

To assess the importance of the bulk non-linear meson interactions on our study on a\ 
and a 2 , we have performed calculations as in Figure 9 but setting r] = 2 with r]i = r\ 2 = 0, 
and setting rji = 1 with r] = 00 and r] 2 = (as indicated, the effect of r] 2 is much smaller 
than that of rji). The results show a completely similar behaviour to Figure 9. Even the 
slope of the contour lines of E s and t in the a\-a 2 plane changes only slightly. Comparing 
with Figure 9, when r] = 2 one finds that E s is shifted by approximately —1 MeV, and that 
when 771 = 1 then E s is shifted by around —3 MeV. The shifts of the surface thickness t are 
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less regular and their magnitude depends on the value of a,\ and d2- 

In order to investigate the impact of the gradient interactions a\ and on the spin- 
orbit potential, in Figure 10 we have plotted V so (z) at the Fermi surface for a few selected 
values of «i and 0:2- The nuclear matter properties and the scalar mass are the same as in 
Figure 6, where we studied the dependence of V so (z) on r] . One can see that the meson 
interaction with «i = 1 and «2 = reduces the strength of V so (z) and shifts the position 
of the minimum slightly to the exterior. On the contrary, the interaction with a\ = 
and a2 = 1 makes the potential well deeper. The combined effect is probed in the case 
a i — a 2 — 1- Since in the relativistic model the spin-orbit force is strongly correlated with 
the Dirac effective mass, we compare in Figure 10 the situation at M^/M = 0.6 and at 
M^/M = 0.7. We realize that the incidence of a\ on V so (z) is weaker for M^/M = 0.7. 
The small perturbations arising from the gradient interactions when M^/M = 0.7 are not 
sufficient to produce a spin-orbit strength equivalent to that of the case M^/M = 0.6. 



3.4 Role of the tensor coupling of the omega meson 

To conclude this section we investigate the influence of the uj tensor coupling 

ifv 



■Pa-VW(z) 



(p a (z) (3.6) 



2M 

which adds some momentum and spin dependence to the interaction. The natural combi- 
nation for this coupling is / v /4. 

Well known from one-boson-exchange potentials (where f v above is commonly written 
as f v /g v ), the tensor coupling was included in the fits to nuclear properties of Refs. P,|T6|j 
(conventional QHD) and ]5||| (effective field theory), and in the study of the nuclear spin- 
orbit force in chiral effective field theories carried out in Ref. |p.7|| . These works noticed the 
existence of a trade-off between the size of the u tensor coupling and the size of the scalar 
field. In other words, the tensor coupling breaks the tight connection existing in relativistic 
models between the value of the nucleon effective mass at saturation and the empirical 
spin-orbit splitting in finite nuclei (which constrains M^/M to lie between 0.58 and 0.64 
M). Including a tensor coupling the authors of Refs. @,f],[n]] were able to obtain natural 
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parameter sets that provide excellent fits to nuclear properties and spin-orbit splittings 
with an equilibrium effective mass remarkably higher (M^/M ~ 0.7) than in models that 
ignore such coupling. We want to analyze the nature of this effect in the simpler but more 
transparent framework of semi-infinite nuclear matter. 

In Figure 11 we have drawn the surface energy coefficient and the surface thickness 
as functions of f v in the range [—0.6, 0.9] for two values of the effective mass and of the 
incompressibility, having set m s = 490 MeV, po = 0.152 fm -3 and a v = —16.42 MeV. To 
exemplify the incidence of f v on the spin-orbit potential, Figure 12 displays V so (z) at the 
Fermi surface for a few of the cases of Figure 11. We also performed the calculations for 
m s = 525 MeV: E s and t are shifted downwards with respect to Figure 11 and V so (z) is 
deeper than in Figure 12, but the global trends with f v are the same. 

Figure 11 shows the strong reduction of E s and t as / v increases (the slope of the 
curves is milder for M^/M = 0.7 than for M^/M = 0.6). Figure 12 reveals that this fact 
is associated with a deeper and wider spin-orbit potential. This agrees with the results of 
Hofer and Stocker p7j who showed in the standard RMF model that the spin-orbit coupling 
reduces the surface energy and thickness. At variance with the individual values of E s and 
t, the ratio E s /t stays to a certain extent constant with f v . 

Figure 12 evinces the sensitivity of V so (z) to f v . The lower the nucleon effective mass is, 
the larger the effect. For M^/M = 0.7 we realize that with positive values of f v (~ 0.3 in 
the present case) one can get a spin-orbit strength comparable, or even stronger, to that 
of the case M^/M = 0.6 and / v = 0, something that could not be achieved with natural 
values of the couplings studied in the previous sections. Since our parametrization with 
M^/M = 0.7 and K = 200 MeV at f v — already has reasonable surface energy and 
thickness (E s = 16.6 MeV and t = 1.97 fm), increasing / v results in smaller values of E s 
and t. This should be compensated with the other couplings (especially a>i and a 2 ) that 
modify the spin-orbit strength to a lesser degree than f v , or the starting point should have 
other values of the incompressibility K and the scalar mass m s . 

The spin-orbit effect has to do with the explicit dependence of the nucleon orbital wave 
functions on the spin orientation A. As described in Ref. |27| nucleons with A = +1 feel 
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an attractive spin-orbit potential and are pushed to the exterior of the surface, whereas 
the spin-orbit force is repulsive for nucleons with A = — 1 which are pushed to the interior. 
As a consequence of this a depletion of particles with A = — 1 occurs at the surface. This 
behaviour is contrasted in Figure 13 for / v = and / v = 0.6 in the case M^/M = 0.7. 
The figure depicts the profiles of the total baryon and tensor densities as well as those of 
their spin components p x (z) and p\{z) for A = ±1, Eqs. ( |2.11| )- (|2.14| ). When the spin-orbit 
strength is large, attraction dominates over repulsion and more particles accumulate at the 
surface than particles are removed from it. Then the total baryon density is enhanced at 
the surface region and it falls down more steeply. 
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4 Surface properties in the asymmetric case 



We briefly recall some basic definitions concerning nuclear surface symmetry properties 
(further details on the relativistic treatment of asymmetric infinite and semi-infinite nuclear 
matter can be found in Refs. p5i^6j^9| ). For a bulk neutron excess 5 = (p n o — Ppo) / (Pno+Ppo) 
(i.e., the asymptotic asymmetry far from the surface), a surface energy coefficient can be 
computed as 

/oo 
dz[£(z)-(a v (5 ) + M)p(z)}, (4.1) 
-oo 

where S(z) is the total energy density of the system of neutrons and protons, a v (5o) denotes 
the energy per particle in nuclear matter of asymmetry 5 , and p(z) = p n {z) + P P {z) with p n 
and p p referring to the neutron and proton densities, respectively. According to the liquid 
droplet model (LDM) for small values of the neutron excess E s (5q) can be expanded 
as follows: 

E M = E *+(%+ 2J t) 62 o + ---- ( 4 - 2 ) 
In this equation J stands for the bulk symmetry energy coefficient, L reads for the LDM 
coefficient that expresses the density dependence of the symmetry energy, and Q is the so- 
called surface stiffness coefficient that measures the resistance of the system against pulling 
the neutron and proton surfaces apart. All of these macroscopic coefficients are familiar 
from semi-empirical LDM mass formulae. 

Another quantity of interest is the neutron skin thickness G, namely the separation 
between the neutron and proton surface locations: 

/oo 
dz[p n (z)/p n0 -p p (z)/ppo] . (4.3) 
-oo 

In finite nuclei would correspond to the difference between the equivalent sharp radii of 
the neutron and proton distributions. In the small asymmetry limit the LDM predicts a 
linear behaviour of G with Sq: 

e = ^*. (4.4) 

For calculations of finite nuclei of small overall asymmetry I = (N — Z)/A, the LDM 
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expansion of the energy can be written as 

'9 J 2 2EX 



E=(a v + JI 2 ) A + 



E K - 



A 2 ' 3 + a c Z 2 A^ 3 + ■ ■ 



(4.5) 



4Q K J 

where ac is the Coulomb energy coefficient. Notice that I ^ 5q in finite nuclei. 

To describe asymmetric matter in the relativistic approach we need to generalize the 
energy density (|2.1| ) by including the isovector p meson. In terms of the mean field R = g p bo, 
with bo the time-like neutral component of the p-meson field, the additional contributions 
to Eq. Q) read 



5>i(*) 



1 



-i4 T3/ 3a.VJ!W 



(p a (z) + -R(z) [p p {z) - p n {z)\ 



2g 2 K K J) 2 



° ' M ) g 2 p 



-Ir 2 (z). 



(4.6) 



The symmetry energy coefficient turns out to be 



J 



+ 



(4.7) 



6 (kl + M*J) 1 / 2 127r 2 m2 1 + t) p {1 - M*JM) ' 

In the conventional model one has f p — rj p — 0. The isovector tensor coupling f p was 
included in the calculations of Refs. p.^|.p|.^6| . The new non-linear coupling rj p between the 
p- and cr-meson fields is of order 3 in the expansion and it has been introduced in Refs. 
Illl . We will not consider higher-order non-linear couplings involving the p meson since the 
expectation value of the p field is typically an order of magnitude smaller than that of the uj 
field [HH . For example, in calculations of the high-density nuclear equation of state, Miiller 
and Serot found the effects of a quartic p meson coupling (i? 4 ) to be only appreciable in 
stars made of pure neutron matter. On the other hand, in analogy to the couplings a\ and 
c*2 for the a and uj fields, we also tested a surface contribution — (Vi?) 2 /(2g 2 ,M) and 
found that the impact it has on the properties we will study in this section is absolutely 
negligible. 

As we have seen, the quantity that governs the surface properties in the regime of low 
asymmetries is the surface stiffness Q. Table 1 analyzes the effect on Q and L of the 
couplings discussed in the preceding sections and of the f p and rj p parameters. On the basis 



of Eq. Q4.4D , we have extracted Q from a linear regression in 5q to fit our results for O up 
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to 5o = 0.1. We have set the equilibrium properties to po = 0.152 fm 3 , a v = —16.42 MeV, 
K = 200 MeV, M*JM = 0.6 and J = 30 MeV, and have used a scalar mass m s = 490 MeV 
and a p-meson mass m p = 763 MeV. Though here we are interested in tendencies rather 
than in absolute values, for comparison we mention that NL1 has (units in MeV) J = 43.5, 
L = 140 and Q = 27 |29|] , the sophisticated droplet-model mass formula FRDM ]34| implies 



J = 33, L = and Q = 29, and the ETFSI-1 mass formula |35] based on microscopic forces 
predicts J = 27, L = -9 and Q = 112. 

Table 1 shows that the influence on the surface stiffness of the volume self-interactions 
?7o, rji and 7/2 is not very large for natural values of these couplings. In the present case Q is 
slightly increased by decreasing f] (i.e., by increasing the quartic vector coupling ( ). For 
7/1 = 1 we find a non-negligible increase of Q, which signals a larger rigidity of the nuclear 
system against the separation of the neutron and proton surfaces. The effect of 772 is again 
moderate as compared to that of 7/1. Q is augmented by a positive r\ p coupling, while a 
negative r] p induces a lower value of Q. Some visible changes in Q take place when the «i 
and a 2 gradient interactions are taken into account. Due to the opposite behaviour of Q 
with cti and ct2, the tendencies compensate in a case like a\ — ol-i — 1, but the net effect is 
reinforced e.g. if «i = — a 2 = 1. 

As one could expect the isoscalar tensor coupling f v has a notable effect on Q, even for 
the relatively small value / v = 0.3 that we have used in Table 1. On the contrary, Q is 
virtually insensitive to the isovector tensor coupling f p . The reason is that the derivative 
of the R{z) field is much smaller than that of the W(z) field. In the least-square fits to 
ground-state properties of Refs. [[§|nj nothing was gained by the p tensor coupling. In any 
case, the best fits of Refs. |5||| have f p ~ 4. 

From Table 1 we recognize that the main changes in the coefficient L arise from the r\ p 
coupling. As a rule of thumb, increasing values of Q are associated with decreasing values 
of L for the bulk couplings 7/ , r/i, 7/ 2 and r\ p . Since L is a bulk quantity, it is not modified 
by the surface interactions. 

Figures 14 and 15 illustrate the dependence on asymmetry of the neutron skin thickness, 
surface energy and surface thickness for some of the cases considered in Table 1. The figures 
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extend up to 5q = 0.3, which widely covers the range relevant for laboratory nuclei (5$ < 0.2). 
As the system becomes neutron rich we can appreciate how a neutron skin develops and G 
grows from its vanishing value at <5 = 0. For small asymmetries the growth is linear in do, 
as predicted by the LDM. The surface energy coefficient grows quadratically with increasing 
neutron excess and the LDM equation (|4.2j) is clearly a good approximation. In general, 
the interactions having thicker neutron skins (smaller values of Q) also have larger surface 
energies. 

The parameter 7] p can be used for the fine tuning of the symmetry properties of the 
interaction without spoiling the predictions for symmetric systems. If in the conventional 
ansatz g p is fixed by the value of the symmetry energy J, in the extended model J depends 
on a combination of g p and tj p , Eq. ( |4.7|) . Therefore, rj p provides in practice a mechanism 
that can help to simultaneously adjust Q (to get the required neutron skin 9) and J (to 
keep the fit to the masses) preserving the symmetric surface properties. 
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5 Summary 



Within relativistic mean field theory, we have investigated the influence on nuclear 
surface properties of the non-linear meson interactions and tensor couplings recently consid- 
ered in the literature. These interactions, beyond standard QHD, are based on effective field 
theories. The effective field theory approach allows one to expand the non-renormalizable 
couplings, which are consistent with the underlying QCD symmetries, using naive dimen- 
sional analysis and the naturalness assumption [|5],|7|)|3|U • 

The quartic vector self-interaction (o makes it possible to obtain a desirable positive 
value of the coupling constant K4 of the quartic scalar self-interaction, for realistic nuclear 
matter properties and within the bounds of naturalness. This (0 coupling has only a slight 
impact on the surface properties. Nevertheless, it helps to find parametrizations where both 
the surface energy coefficient E s ant the surface thickness t lie in the empirical region. The 
Co vector non-linearity makes the spin-orbit potential well deeper, although the effect is 
almost negligible. Concerning the volume non-linear couplings rji and 772, they also allow 
one to obtain positive values of K4 in the region of naturalness, depending somewhat on the 
saturation properties (incompressibility and effective mass). The surface properties are not 
much affected by these bulk terms either, and it turns out that 772 has a marginal effect as 
compared to that of 771. 

The equilibrium properties do not depend on the couplings a% and that involve the 
gradients of the fields. Thus, these couplings serve to improve the quality of the surface 
properties without changing the bulk matter. In the conventional a — u model the only 
parameter not fixed by the saturation conditions is the mass of the scalar meson. In the 
cn-\—Oi2 plane the lines of constant E s have a different slope than those of constant t. It is 
then possible to keep a fixed value of E s and to modify the value of t by choosing a\ and 
a?2 appropriately. The range of variation of E s and t with ot\ and oli is wider than with the 
volume couplings. This justifies including these gradient terms in the energy functional in 
spite of being of order 5 in the expansion. The aii and ai surface meson interactions also 
influence the spin-orbit potential, but the effect is not extremely significant. 

The effective model is augmented with a tensor coupling of the u meson to the nucleon. 
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An outstanding feature is the drastic consequences it has for the spin-orbit force. We have 
emphasized how inclusion of f v permits to obtain a spin-orbit strentgh similar to that of 
M^/M ~ 0.6 with larger values of the equilibrium nucleon effective mass, contrary to the 
phenomenology known from models without such a coupling. 

We have discussed the implications of the extra couplings of the extended model on 
various surface symmetry properties. We have restricted ourselves to the regime of low 
asymmetries, where the liquid droplet model can be applied and the surface stiffness coeffi- 
cient Q is the key quantity. In particular we have pointed out the role that the non-linearity 
i] p of the isovector p-meson field may play in the details of the symmetry properties. 
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Table captions 

Table 1. The surface stiffness coefficient Q and the coefficient L for several values of the 
couplings analyzed in the text. We have set p — 0.152 fm~ 3 , a v = —16.42 MeV, 
K = 200 MeV, M*JM = 0.6, J = 30 MeV and m s = 490 MeV. 
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Figure captions 



Figure 1. The couplings ^ s /4vr, g v /4:ii, k 3 , k 4 and Co versus the parameter 770 defined in Eq. 
(|3.3| ). The naturalness assumption requires all these couplings to be of order unity. 
We have taken p = 0.152 fm~ 3 (k F = 1.31 fm" 1 ), a v = -16.42 MeV and m s = 490 
MeV. 

Figure 2. Surface energy coefficient E s for several values of the parameter rj , K, M^/M 
and m s , with p = 0.152 fm -3 and a v = —16.42 MeV. 

Figure 3. Surface thickness t of the baryon density profile for several values of the param- 
eter 770, K, M^/M and m s , with p = 0.152 fm" 3 and a v = -16.42 MeV. 

Figure 4. Surface energy coefficient E s and surface thickness t for several values of the 
parameter r/ , M^/M and m s , with p = 0.1484 fm" 3 , a v = -15.75 MeV and K = 200 
MeV (Ref. [|H|). 

Figure 5. Baryon density p(z) and surface tension density a(z) = S(z) — (a v + M)p(z) 
of semi-infinite nuclear matter for some values of the parameter 770. It is po — 0.152 
fm" 3 , a v = -16.42 MeV, K = 200 MeV and m s = 490 MeV. 

Figure 6. Orbital part of the spin-orbit potential V so (z) at the Fermi surface and central 
mean field V c (z), Eqs. (|2.17|) and (|2.19|) respectively, for some values of the parameter 
r] . It is po = 0.152 fm" 3 , a v = -16.42 MeV, K = 200 MeV and m s = 490 MeV. 

Figure 7. The couplings g s /47r, g v /4ir, k 3 and K4 against the parameters 771 (left) and 772 
(right). With p = 0.152 fm" 3 , a v = -16.42 MeV and m s = 490 MeV, results are 
shown for K = 200 MeV and M^/M = 0.6 in part (a), and for K = 350 MeV and 
M^/M = 0.7 in part (b). 

Figure 8. Level curves in the plane r/i-r/2 of the surface energy coefficient E s (in MeV, solid 
lines) and of the surface thickness t (in fm, dashed lines), with ( = a\ = a 2 = f v = 0. 
The point 7/1 = r/2 = is marked by a cross. It is po = 0.152 fm -3 , a v = —16.42 MeV, 
K = 200 MeV, M^/M = 0.6 and m s = 490 MeV. 
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Figure 9. Same as Figure 8 in the plane ai~a 2 , with Co = Vi — V2 = fv = 0. The point 
&i = ol 2 = is marked by a cross. 

Figure 10. Orbital part of the spin-orbit potential V so (z) at the Fermi surface for some 
values of the couplings ai and a 2 - The equilibrium properties of nuclear matter and 
the scalar mass are the same of Figure 6. 

Figure 11. Surface energy coefficient E s and surface thickness t as functions of the strength 
/ v of the cu-meson tensor coupling, with Co = ^1 = ^2 = «i = «2 = 0. We have set 
p = 0.152 fm- 3 , a v = -16.42 MeV and m s = 490 MeV. 

Figure 12. Orbital part of the spin-orbit potential V so (z) at the Fermi surface for some 
values of the tensor coupling f v . The equilibrium properties of nuclear matter and the 
scalar mass are the same of Figures 6 and 10. 

Figure 13. Total baryon density p(z), total tensor density Pt(z) and their components 
p x (z) and Pt(z) for the two spin orientations A = ±1. They have been calculated for 
f v = and / v = 0.6, with M*JM = 0.7, K = 200 MeV and m s = 490 MeV. 

Figure 14. Neutron skin thickness O as a function of the bulk neutron excess 5 - The solid 
line is the result of the conventional model (Co = Vi = V2 = Vp = a i = a 2 = fv = f P = 
0). The other lines differ from the latter in the indicated parameter. 

Figure 15. Same as Figure 14 for the surface energy coefficient E s and the surface thickness 
t as functions of the bulk neutron excess squared 5q. 
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